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I ■ Abstract 

C . In order to clarify how the percolation theory governs the conductivities in 

Q \ real materials which consist of small conductive particles, e.g.^ nanoparticles, 

with random configurations, we numerically investigate the conductivities of 
CN ■ continuum percolation models consisting of overlapped particles using the 

^ ■ finite difference method as a sequel of our previous article (Int. J. Mod. 

\/^\ Phys. 21 (2010), 709). As the previous article showed the shape effect of 

each particle by handling different aspect ratios of spheroids, in this arti- 
cle we numerically reveal influences of the agglomeration of the particles on 
Q ■ conductivities after we model the agglomerated configuration by employing 

a simple numerical algorithm which simulate an agglomerated configuration 
of particles by a natural parameter. We conclude that the dominant agglom- 
eration effect on the conductivities can be interpreted as the size effect of an 
/\ \ analyzed region, and shape of the agglomerated clusters also has an effect on 

c^ ■ its universal property. 

Key words: 

continuum percolation, conductivity, critical exponents, agglomerated 



^Present address: Microcreate System Co., 2-1-15 4C, Chuou Yamato Kanagawa 242- 
0021 Japan. 



Preprint submitted to arXiv March 16, 2012 



cluster, conductive nanoparticle, finite difference method, 



1. Introduction 

As in the previous article [l| , we investigate conductivities of a percola- 
tion model using the finite difference method (FDM) |2|. The purpose of this 
article and the previous one is to clarify how the percolation theory governs 
the conductivities in real materials which consist of small conductive parti- 
cles, e.g., nanoparticles, with random configurations. In general, the small 
particles in a real material have some typical sizes, shapes like spheroids, 
and cluster structures due to their agglomeration (3|-J8|. The purpose of this 
article is to investigate an agglomeration effect on the conductivity whereas 
the previous one [l| is a study on the effect of the shape of each particle. We 
propose a simple numerical algorithm to simulate an agglomerated configu- 
ration of particles though there are some phy sical models which numerically 
simulate the agglomerated phenomena |9|, [lO|. In this article, we employ the 
simple algorithm. Since our algorithm has a parameter which controls the 
agglomerated configurations consistently without any difficulties, we numeri- 
cally show how the conductivity depends upon the agglomerations. Then we 
reveal the fact that the main agglomeration effect on the conductivity can be 
interpreted as the size effect of the system and the shape of the agglomerated 
cluster also has an effect on the universal property of the conductivity. 

The conductivity is one of the most concerned phenomena in the per- 



colation theory [ll[-[l6| and it is important to investigate the conductivity 
curve 

where c is a constant factor, p is the volume fraction, pc is the percolation 
threshold, and t is the critical exponent of the conductivity. The behaviors 
of the threshold pc and the exponent t represent the conductive properties 
of the systems in the percolation phenomena, which have some universal 
properties. It means that their dependence on some parameters of a model is 
sufficiently weak due to the randomness. In order to discuss the universality, 
the percolation theory is, ordinarily, based upon a system with the infinite 
size 
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However when we apply the percolation theory to the real composite 
materials of conductive nanoparticles, we encounter effects coming from the 
shape of the nanoparticles and several characteristics lengths. Besides lattice 



percolation models which are given on discrete lattices, the continuum per- 



colation model (CPM) was introduced as in Refs. [ll|, p. 108-111] and [17|. In 
order to handle the size and the shape of particles, CPM has been studied. 
In the previous article [l| , we studied the shape effects on the conductivities 
in CPMs among the different aspect ratios of spheroids, where we allowed 
overlap of the particles. The previous article shows that the conductivity in 
CPM strongly depends upon the shape of the composed particles. 

In this article, we will consider another shape effect or the agglomer- 
ated effect of these stuffed spherical particles on the conductivity. With the 
development of the technology, the smaller the size of the (nano-)particles 
becomes, the smaller the size of the related devices becomes. When we in- 
vestigate the conductivities in the real composite materials of conductive 
nanoparticles from the viewpoint of the percolation theory, we handle three 
characteristic lengths, i.e., 1) the size of the particle as a minimal size, 2) 
the size of the system, e.g., the thickness of the film, as the maximal size, 
and 3) the size of the percolation cluster which is given by a multiplication 
of the size of the particles and is also related to the size of the system. When 
the size of devices is sufficiently small, it is important to consider the size 
effect of the system though the conductivities in the system with infinite size 
are basically concerned in the percolation theory. In other words, in order 
to apply the percolation theory to a real material system, it is crucial to 
consider the relations among these scales. 

The smaller the size of particles is, the more agglomerated the particles 
becomes, due to their interface energy. Agglomeration forms agglomerated 
clusters and the cluster brings the fourth characteristic scale to the system. 
Thus the evaluation of the agglomeration effect on the conductivity is very 
important. By employing the simple algorithm to simulate the agglomeration 
of particles, we numerically give a series of agglomerated configurations of 
CPM, which we call agglomerated continuum percolation models (ACPMs), 
and investigate the conductive properties of ACPMs by means of FDM. In 
this article, we also handle the overlapping particles [l| . We show the depen- 
dence of the conductivity curves on the agglomeration as in Figure [71 

Since the finite size effects in the conductivity on a percolation model 



were studied in Chapters 4 and 5 of Ref. [11|, based upon these studies, 
we numerically consider the relation between the geometrical effects and the 
properties of the conductivities in ACPMs. As a result, we show that one of 
the dominant agglomerated effects on the conductivities with ACPM can be 
regarded as the size effects of the system and the shape of the agglomerated 



clusters also affects the conductivity. The shape effect implies that the ag- 
glomeration would have an effect on the universality class of the conductivity 
curve ([1]) , and thus we numerically show the fact in Figure [T2j 

Contents in this article are as follows. Section 2 shows our computational 
method. Subsection 2.1 describes our algorithm to construct the agglomer- 
ated clusters in CPM and the geometrical setting in ACPMs. Section 2.2 
provides the computational method of the conductivity over ACPMs using 
FDM, which is basically the same as that in the previous article [jj . Section 
3 shows our computational results of the conductivity in ACPMs with the 
agglomerated clusters. In Section 4, we discuss our results from geometrical 
viewpoints. In Section 5, we summarize the results and the discussions. 

2. Geometrical setting of ACPM 

In this section, we show our geometrical setting of ACPM. We model 
the agglomerated clusters out of a simple algorithm which is governed by a 
parameter 7agg G [0, 1]. We also briefly show the computational method of 
the conductivities over ACPMs usingFDM in Subsection 2.2, whose details 
are described in the previous article [l|. 

2.1. Agglomeration algorithm 

We set particles parametrized by their positions {x, y, z) into a box-region 
B := [0, Xo] X [0, Ho] X [0, 2:0] at random and get a configuration 9^„ as one of 
CPMs. In this article, we set Xq = z/q = Zq = L. The particle corresponds 
to a stuffed sphere or ball with the same radius p, B^^^y^^z, '■= {{x,y,z) e 
B I \{x,y,z) — {xi,yi,Zi)\ < p}. The configuration £H„ is given by 9^„ := 

By fixing the radius p of the particle, p = 1, and a number 7agg G [0, 1] 
which is called agglomeration parameter, we introduce an algorithm to con- 
struct the configuration 9^„ in ACPM. 

We illustrate our algorithm by a flowchart in Figure [1] As an initial state, 
the configuration 9^o has no particle. As the first step, for a uniform random 
position (x, y,z) E B we set a particle B^y^^ whose center is (x, y, z) and the 
radius is p, i.e., 9^i := B^^y^- 

Let us consider n-step. We take a position (x, y, z) at uniform random 
in B, and another random parameter 7 at uniform random in [0, 1]. If the 
random parameter 7 is greater than 7agg, we employ the position (x, y, z) as 




Stiirt 



Set a ball at a random position (x.j',z) 



Pick up a random number /. 



n ^n+\ 




Get a. random position (a-j;s) . 



SeL a ball at (x^v.z) 




Get a random positiori (xj\z) . 




Yes 



No 



Figure 1: The flowchart of the agglomeration configuration algorithm. 



the center of a particle and 9^n+i := ^U-^a;,j/,2- ^^ i^ noted that we allow 
the particles to overlap each other. 

When the random parameter 7 is given as 7 < 7agg, we first check whether 
the ball whose center is the position [x, y, z) is connected with the previous 
configuration 9^„ or not. If it is connected with the configuration £H„, or 
^n n ^x,y,z 7^ 0, we employ the position and add the particle into the con- 
figuration 9^„, or 9^„+i := 1H„ |J-Sx,y,^- 

If not or 9^„ f] -Bx,j/,z = 0, we abandon the position and go on to take 
another uniformly random position (x, ?/, z) in B until we find the position 
which supplies a connected particle B^^y^z with 9^„. 

In other words, if we take 7 which is smaller than 7agg, the added particle 
must be connected with the previous configuration 9l„. Thus, 7agg stands for 
the agglomeration of the particle system. 

By monitoring the total volume fraction which is a function of £H„ and is 
denoted by vol(£H„), we continue to put the particles as long as vol(9^„) < p 
for a given volume fraction p. We find the step n{j)) such that vol(9^.„(p)_i) < 
p and vol(£H„(p)) > p. Since the difference between \o\{^n(p)-i) and vol(9^„(p)) 
is at most 0.9 x 10~^ for the employed parameters. Thus we regard vol(£H„(p)) 
as the volume fraction p itself hereafter under this accuracy. 

As illustrated in Figure 121 the uniform random configuration is realized 
as the case 7agg = 0, and the agglomeration configurations are given with our 
agglomeration parameter 7agg > 0. We regard ACPM of the case 7agg = as 



the ordinary CPM for balls with the same radius [ll|, ll7| and thus we simply 
call it CPM hereafter. 

Since we use the pseudo-randomness to simulate the random configuration 
9i„(p) for given p and 7agg, the configuration ^n{p) depends upon the seed ig of 
the pseudo-randomness which we choose. We let it be denoted by ^7^gg,p,j^. 
For the same seed is of the pseudo-randomness, a configuration 9^^^ p j^ of a 
volume fraction p naturally contains a configuration 91y^gg,p',i^ for p' < p due 
to our algorithm, i.e., 9^7^ggy,j^ C 9iy^gg,p,j^. Hence the elements in the set of 
the configurations {9^7^gg,p,J3 | p G [0, 1]} with the same seed is are relevant. 
As shown later, the conductivity curve ([T]) over {D^-y^gg^p^j^ \ V ^ [0)1]} fo'^ 
fixing is is obtained as a smooth curve even for the non- vanishing 7agg- 

In order to illustrate this algorithm, we will show two dimensional case. 
In a vary similar way to the three-dimensional case, we have two dimensional 
agglomeration configurations which are given in Figure |3l Figure H] exhibits 
the transmission electron micrograph (TEM) of the real agglomerated clus- 



ters in Refs. [8| and [20| . Ref . |8| is of the composite of polyacrylonitrile with 



Volume Fraction p=0.2 




(a) 








(h) 







(c) 



(d) 



Figure 2: The agglomeration parameters: Random particle systems with the agglomera- 
tion parameter 7agg = 0.0,0.5,0.8 and 0.95 for (a), (b), (c) and (d) respectively. 



acetylene black particles, and Ref. [20| shows the epoxy nanocomposites with 
silicate oxide particles. These configurations in Figure E] might simulate some 
agglomeration states in Figure HI Even though it is difficult to visualize the 
three dimensional cases well, it is expected that our algorithm generates the 
similar agglomerate configurations in Figure |2j 



Volume fraction p = 0.4 




Figure 3: The agglomeration parameters: Random particle system with the agglomeration 
parameter 7agg — 0.0,0.5,0.8 and 0.95 for (a), (b), (c) and (d) respeetively. 



2.2. Computation of the conductivity in ACPMs 

To apply FDM [2| to the computation of the conductivity in ACPM, 
we use a A''^. x Ny x N^ lattice denoted by C to represent the box-region 
B by the integers A*"^, Ny and Nz- In this article, we mainly assume that 
Nx = Ny = Nz = 216 and the radius of the particle p = 1 corresponds to six 
meshes. It means that the ratio between the volumes of B and a particle is 
about 1.1 X 10^ and the size oi B = [0, L]^ is given as L ^ 36. 

Further we set the conductivity distribution o"(x, y, z) which consists of 
the conductive particles and the insulator as a background in the box-region 
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(a) 



(b) 



Figure 4: TEM micrograms of (a) polyacrylonitrile nanocomposite of acetylene black 
particles in Ref. [8| and (b) epoxy nanocomposites containing silicate oxide particles in 
Ref. ii. 



B as in Figure [21 We put the conductivity density (Tmat = 1 inside of each 
5^7agg,p,is- As the insulator, we set the infinitesimal conductivity Uinf = 10~^ 
within background. In other words, we handle binary materials with largely 
different local conductivities o"mat and ainf. 

In order to compute the total conductivity, we set the voltage = 0o = 1 
and = on the upper and the lower faces, i.e., [0, Xq] x [0,yo] x {zq} and 
[0, Xq] X [0, yo] X {0} respectively as the boundary condition corresponding to 
the electrodes. As the side boundary condition, we used the natural boundary 
for each yz-ia.ce and each xz-iace. In other words, at the side boundaries, we 
imposed that current normal to each face vanished. 

Following our algorithm of FDM as mentioned in detail in Ref. [l|, we 
numerically solved the generalized Laplace equations over the region, 

V ■ (tV0 = 0, (2) 

for the conductivity distribution a{x,y,z) for each 9^^^ pj^ to obtain the 
potential distribution (j){x,y,z). By numerically solving (|2]), we obtained 
the total conductivity Utotai of the system after we integrated the current 
(jV0 over a xy-plane. Since cxtotai is determined for each 9^-y^gg,p,j^, crtotai 
is a function of the volume fraction p, 7agg and the seed ig of the pseudo- 
randomness. We denote it by crtotai(p,7agg, «s) explicitly or simply (Ttotai(p)- 

As mentioned in Introduction, the conductivity curve crtotai(p) {p G [0, 1]) 
is expressed by 



(p-PcY 



ioT p e [pcA], 



crtotM = { (^•°-^^)*' ^-;/^-.L^-^J' (3) 

otherwise, 

where Pc is the threshold and t is the critical exponent, or merely exponent. 
Since the threshold pc and the exponent t are determined by {9^7agg,p,is I 
p G [0,1]}, they are functions of the agglomeration parameter 7agg and 
the random seed is- Thus we sometimes express them as pdlaggy'is) and 

We evaluated the threshold Pc(7agg5 "^s) and the exponent t{pfa.gg^ is) as the 
fitting parameters so that each average of the square error from the curve 
is the smallest, by monitoring the square root of the average of the square 
error. The square root of the average of the square error is denoted by (5crtotai 

= <^0-total(7agg,'^s)- 
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3. Results 



As mentioned above, for each seed ig of the pseudo-randomness, a config- 



uration IH, 



,p,^i 



to a volume fraction p naturally contains the configuration 



91. 



,p ,«£ 



to p' < p of the same seed is- Hence for a given seed is, the to- 



tal conductivities atotai between 9^7^gg,p,j^ and D\. 
we obtained the conductivity curve over {Dl. 
curve even with the non- vanishing 7agg. 
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i^ are relevant. Thus 
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Figure 5: The conductivity curves for agglomeration parameters 7agg = 0.0, 0.5, 0.8 and 
0.95 of the random seed ig — 1: (a) is of the Unear scale and (b) is of the logarithm scale. 

Figure [5] illustrates the conductivity curves for the seed is = 1- Figure 
Mjo) shows that we handled the binary conductive materials. When p = 0, 
the total conductivity Utotai is equal to Uinf which stands for the insulator 
and is negligible if we dealt with atotai linearly as shown in Figure [5](a). 
By using the least square error fitting method, we evaluated the threshold 
Pcilagg, is) and the exponent t(7agg, is) from the curve. In the evaluation, we 
used the linear scale of atotai or the curves in Figure [5](a). Figure [6] shows the 
fitting errors (5o"totai v.s. the seed ig] 5o"totai is the average of the squares of the 
deviations from the curve ([3]) over the fitting points for each 7agg and is] the 
fitting points are 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.6, 0.7, 0.8, and 
0.9. We computed thirty curves for each case, 7agg = 0.0, 0.5, 0.8 and 0.95. 
The deviations Satotai are not large even for non-vanishing 7agg compared 
with o"mat = 1- The distribution of ^Utotai shows that the larger 7agg is, the 
larger the values are. Its maximum is 5crtotai,max < 0.0084 for 7agg = 0.95, but 
the average ^cxtotai is less than 0.0048; the values of (7agg, ^o'totai) are given as 
(0.0,0.0035), (0.5,0.0043), (0.8,0.0047) and (0.95,0.0048). 



11 





0.012 




0.01 




0.008 


B 
o 


0.006 




0.004 




0.002 








Y,gg = 0.50 --X- 

Tagg = 0-80 *■ 

0.95 B 




10 



15 

Seed L 



20 



25 



30 



Figure 6: The fitting errors (5crtotai vs the random seed is for 7agg — 0.0, 0.5, 0.8 and 0.95. 



It means that the thresholds Pc and the exponents t represent the con- 
ductivity curves well in our computations. In fact since we have 



9crtotal(p) 



-O-total(p)- 



1 — p 



datote,i{p) , ., P-Pc 

= CTtotalbJlog- 



the accuracies 6aP and 6at of pc and t, could be estimated using ^crtotai; 



5aP 



da 



total 



dpc 



6a, 



total ! 



SJ 



da 



total 



dt 



6a, 



total • 



For example, around p ~ 0.5 (atotai ~ 0.2), they are evaluated as S^p ~ 
(Jatotai ~ 0.005 and 5^t ^ 4.35owi ~ 0.022 by setting pc ~ 0.27, t ~ 1.63 as 
shown below. 

As the case 7agg = 0, our computational results are pc = 0.273 ± 0.012 
and t = 1.628 ± 0.036. It should be noted that these computational results 
are obtained by FDM method with finite lattice. As the previous article, we 
showed a finite size effect of the lattice by using an extrapolation scheme of a 
lattice-size dependence of pc and t, these computational results pc = 0.273 ± 
0.012 and t = 1.628 ± 0.036 are on the extrapolation line in Ref. |l|; The 
extrapolated values as shown in the previous article [1| are pf^ = 0.296±0.013 
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and t^^ 
and t 



1.580±0.042, which agree with p^ = 0.289573±0.000002 in Ref. [23 



1.6 ± 0.1 in Refs. 21 



and 



22l . Hence our computational scheme is 



m 23 



consistent with these results l|, |21 

In this article, since we focus on the difference of the conductive properties 
among several 7agg's, we basically fix the lattice-size and investigate the con- 
ductivities in the finite region B without any corrections for the finite-lattice 
effects. 
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(a) (b) 

The agglomeration dependence of the thresholds (a) and critical exponents (b). 
circle corresponds to each computation for a seed and the dotted line shows the 



The dependencies of exponents and thresholds upon the agglomeration 
parameter 7agg are illustrated in Figure [7] and Table [1] Figure [7] shows that 
the agglomeration parameter 7agg has large effects on the conductivity over 
ACPM beyond the accuracies 6„Pc ~ 0.005 and 6„t ~ 0.022. The larger 
the agglomeration parameter 7agg is, the more largely the threshold and the 
exponent depend on the seed is- It means that the fluctuation coming from 
the randomness is enhanced by the agglomeration parameter 7agg- 

Furthermore the larger the agglomeration parameter 7agg is, the smaller 
the trend of the threshold is and the larger that of the exponent is. As shown 
in Table [H the larger the agglomeration parameter 7agg becomes, the smaller 
the average of the thresholds p^(7agg) is and the larger the average of the 
exponent t(7agg) is. 
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Tabic 1: The 7agg dependence of the threshold and the exponent. 





Threshold 


Exponent 


7agg 


Average 


Maximum 


Minimum 


Average 


Maximum 


Minimum 


0.0 


0.273 


0.284 


0.261 


1.628 


1.661 


1.588 


0.5 


0.250 


0.283 


0.231 


1.689 


1.758 


1.597 


0.8 


0.224 


0.324 


0.159 


1.772 


1.951 


1.506 


0.95 


0.202 


0.276 


0.047 


1.809 


2.214 


1.523 



4. Discussion on the agglomeration effects 

We now consider the reason why the agglomeration has the effects on the 
conductivity in our ACPM. 

We have dealt with the percolation phenomena on a finite region B 
whereas it is well-known that the conductivity has dependence on the size of 



region as in Chapter 5.1 in Ref. [ll|. We naturally have the characteristics 



length Tc induced from the size of the particles or the radius of the particle 
p; for vanishing 7agg, i.e., Tc(7agg = 0) = p. Due to the finiteness of size 
L of B, the ratio L/vc has an effect on the conductivity; at 7agg = 0, we 
have that L/vc ~ 36. For non- vanishing agglomeration parameter 7agg 7^ 0, 
it is expected that the characteristic length Tc is larger than that of CPMs 
(7agg = 0) due to the agglomeration, i.e., Tc(7agg 7^ 0) > p. Further by let- 
ting ^ = ^(p, 7agg) be the correlation length which represents the percolation 
phenomenon, we should compare the size of region L, the correlation length 
^ and the characteristic length re 
total conductivity behaves like 



According to Chapter 5.1 in Ref. [ll|, the 



0-total(P, 7agg) OC 









for L/vc > (^/rc 
for L/vc < (^/rc 



(4) 



where p is the non- negative parameter (73a) in Ref. [ll[. Since it is known 
that ^ diverges at the critical point p^ and in every numerical estimation 
we handle only a finite region, the formula dl]) means that every numerical 
estimation gives higher conductivity in the vicinity of the point Pc than that 
of an infinite region apart from the random fluctuation. 

Though the percolation theory is basically concerned with conductivity 
curve in the infinite region, we concern ourselves about the size effect of the 
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real materials. Further in Ref. [1|, we have showed that the shape effect 
has crucial effects on the conductivity curves. Thus we will consider the 
geometrical properties of agglomerated clusters. 

4.I. The geometrical features in ACPM 

Let us consider the size and the shape of the agglomerated clusters which 
depend on the agglomeration parameter 7agg. 

4-1.1. The characteristic length r^. in ACPM at p = 0.1 

First we consider the size effect and give an estimation on the character- 
istic length Tc for non- vanishing 7agg. Since it is difficult to estimate it, we 
compute the difference among the size of the isolated agglomerated clusters 
at a lower volume fraction p than pc- In other words, we compute statistical 
behavior of the percolation clusters, or the connected particles, at p = 0.1. 

When the agglomeration parameter 7agg becomes large, it is expected that 
the size of the percolation cluster becomes larger than one of the uniform 
randomness or the case (7agg = 0). The size of the percolation cluster is 
associated with the correlation length ^(p, 7agg) in (^. 

In order to evaluate the effect of the size, we consider effective radius 
Pcist and maximum length Lmax of each agglomerated clusters in ACPM, and 
the numbers A^'agg of agglomerated clusters. Here the effective radius pcist is 
defined such that inpli^^/S is equal to the volume of the percolation cluster. 
First, we consider a histogram of the effective radius Pcist- Figure |H] shows 
frequency of the effective radius pcist at the volume fraction p = 0.1 which is 
smaller than any thresholds p^(7agg)- Figure [8] means that the larger the 7agg 
is, the larger the size of the cluster is. 

Figure |8] shows the fact that the histogram has two peaks; the first peak 
at Pcist = p as the size of isolated particle, and the second peak around 
Pcist = v^ ~ 1.26p as the size of two particles. The probability of two 
slightly connected particles is larger than the probability of the state which 
has Pcist £ (p, v^p) because of radial measure Pcist^^Pcist for the radius pdst 
from the center of a particle. The percolation is determined by the largest 
cluster size but it should be statistically treated. We consider the right hand 
side of the second peak of the histograms. We fit the shape of the histogram 
over Pcist G [1-3,8] by /i(pcist) = Aexp(-(pcist - p)/Pcist) well using the least 
mean square error method, where A and p^i^^ are fitting parameters. Then 
p°ig^ is given in Table [21 Here 6h is the square root of the average of the least 
mean square error. 
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Figure 8: The histogram of size pcist of percolation cluster for each agglomeration param- 
eter 7agg at p = 0.1. 



Table 2: The pdst dependence on the agglomeration parameter 7, 



7agg 





0.5 


0.8 


0.95 


pUp 


0.315 


0.458 


0.476 


0.522 


6h 


0.020 


0.019 


0.024 


0.028 
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Table [2] shows that the agglomeration parameter 7agg increases the size 
of the clusters at p = 0.1. The characteristic length r^ is directly relevant 
to the size of clusters p^isf Though it is difficult to evaluate the difference 
of the size of clusters p^j^^ and also Tc at p > pc, it is expected that it plays 
similar roles even for every p G [0, 1]. 



Table 3: The Pcist, ^max, and A'agg of the agglomerated (pcist > 1-3) particles dependence 



on the agglomeration parameter "fagg at p = 0.1. 



7agg 





0.5 


0.8 


0.95 


Pclst 


1.572 


1.800 


1.978 


2.200 


-^max/ ^ 


2.803 


3.310 


3.668 


4.009 


iVagg 


2977 


2725 


1722 


952 



Further Table [3] shows the averages of the size pdst and the maximum 
distance L^ax of the percolation clusters which have pcist > l-3p, and their 
number iVagg. The larger 7agg is, the larger Lmax and p^st are, and the smaller 
A^agg is. 

They mean that the larger 7agg is, the smaller the ratio L/vc is. Hence it 
is expected that the larger 7agg is, the larger the fluctuation becomes. 

4-1.2. The shape of the agglomerated cluster in ACPM 

Here we consider the shape of the agglomerated cluster. In our algorithm 
which is shown in Figure [H a new particle with 7 < 7agg must be connected 
with the particles which have been already placed. If the volume fraction 
is much less than the threshold, the connected (agglomerated) cluster which 
consists of A^ particles could be regarded as an orbit of a random walk for 
discrete N time step. 

In fact, the size of the agglomerated cluster is proportional to the square 
root of A^ as shown in Figure [9l Figure [9] displays the correlation between the 
Lmax and yN for 7agg = 0.0, 0.5. They are linearly relative. This coincides 
with the properties of the random walk due to the central limit theorem. 

Figure [9] means that if we regard the agglomerated cluster as a cylinder, 
the radius around the center axis is proportional to V^ since the volume 
should be proportional to A^ but Lmax oc a/ZV. For a sufficiently large A^, 
the agglomerated cluster might be regarded as a thin cylinder rather than a 
thick cylinder. 
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Figure 9: The dependence of the Lmax on the ^/N for 7agg = 0.0, 0.5. 



4-2. The agglomeration effect on conductivity in ACPM 

In this section, we will investigate the behavior in Figure [7] as the agglom- 
eration effect on conductivity in ACPM. 

As shown above, when we fix B, the larger rc(7agg) becomes, the smaller 
the effective size L/vc in @ is regarded. It is expected that the total conduc- 
tivity fTtotai(p,7agg,^s) strongly depends on the configuration ^j,^g,p,i, since 
statistical average does depend upon the ratio L/r^ which is not sufficiently 
large for 7agg > 0. As shown in Figure [3 it means that the larger 7agg is, the 
more largely the threshold Pc(7agg 7^ 0,2^) and the exponent t(7agg 7^ 0,is) 
depend on the seed ig. Thus we consider the size effect first. 

4-2.1. The size effect on conductivity in CPM 

In this subsection we consider the size effect on conductivity in CPM or 
the case 7agg = by changing the radius p directly. 

We computed the threshold Pc and the exponent t for the different radius 
p(> 1) by fixing the size of B. (From the numerical viewpoint, we performed 
the similar computations of 7agg = for small N^ = Ny = N^ by fixing the 
mesh of p. It corresponds to the variation of the radius p of the particles for 
fixing B relatively.) 

Figure [TU] and Table H] give the dependence of the threshold and the 
exponent on the size of particles for the same B. Figure [10] shows that the 
(relatively) larger the radius p is, the larger the fluctuations of the threshold 
Pc and the exponent t are. This property is similar to Figure [3 However 
the fluctuations of both threshold and exponent in Figure dU] look symmetric 
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Figure 10: Thresholds (a) and critical exponents (b) vs radius of particles in the same B. 
The fiUed circle corresponds to each computation for a seed and the dotted line is their 
average. 

with respect to their averages whereas the fluctuations in Figure [7] show 
asymmetry. 
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Figure 11: Critical exponents vs thresholds (a) for various p with 7agg 
various 7agg with p = 1. 



0.0, and (b) for 



The equation (j4]) is evaluated from the statistical viewpoint. There is no 
difference among p's if B has inflnite region. Since B has a flnite size, for 
large p, we don't have so sufficiently many particles in B that its statistical 
average works well. Hence the fluctuation is enhanced. In other words, the 
dependence of the conductivity curve ([3]) on the random seed is is larger than 
that of p = 1 as in Figure [TDl 
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Since (jlj) means that the finiteness of B makes the threshold p small sta- 
tistically. In other words, the trend of the statistical average of the threshold 
shows that the larger the size p in CPM is, the smaller the threshold %{p) is 
as shown in Table [Hand Figure [TO] (a). Here, %{p) stands for the statistical 
average of the pdp, is) over the random seed ig. This property of the thresh- 
old is also similar to Figure [7] (a), though the dependence of the exponents 
in Figure [TUJ (b) is quite different from Figure [7] (b). 

When the individual total conductivity <7tota.\{p, P, is) has smaller pc{p,is) 
than %{p), the non-vanishing crtotai(p,p,«s) at p € iPc{p,is),%{p)), must 
increase weakly with respect to p. It means that the exponent t{p, is) becomes 
larger in the conductivity curve (^ than t. It implies that the exponent 
t{p, is) and the threshold Pc{p, is) in CPMs of different p and is are correlative 
and that the smaller the thresholds Pc{p, is) are, the larger the exponents 
t{p,is) are. Particularly Figure [TTT a) exhibits the correlation between the 
exponent t and the threshold Pc in CPMs. 

Table 4: The p dependence of the threshold and the exponent. 





Threshold 


Exponent 


p 


Average 


Maximum 


Minimum 


Average 


Maximum 


Minimum 


1.0 


0.273 


0.284 


0.261 


1.628 


1.661 


1.588 


1.35 


0.269 


0.291 


0.246 


1.620 


1.711 


1.530 


2.0 


0.274 


0.305 


0.232 


1.572 


1.742 


1.457 


2.7 


0.260 


0.324 


0.192 


1.581 


1.836 


1.340 


3.86 


0.242 


0.345 


0.152 


1.568 


1.873 


1.253 



4-2.2. The size effect on conductivity in ACPM 

Following the above discussions, we consider the agglomeration effect on 
the conductivity in our ACPM. 

Subsection 14.11 means that the larger the agglomeration parameter 7agg 
is, the larger the size of the percolation clusters and characteristic length Tc 
are and the smaller the effective size L/tc in (jl]) is. For the agglomeration 
parameter 7agg deviated from 0, B corresponds to a relatively smaller region 
than the uniform random case (7agg = 0). Subsection 14 . 2 . 1 1 means that due to 
the finite size effect, the dependence of the threshold and the exponent upon 
the seed of the pseudo-randomness, i.e., these fluctuations, are larger than 
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the case of 7agg = 0. In other words, the large fluctuation for non- vanishing 
7agg comes from the finiteness of B and the size of the agglomerated clusters. 
Since the agglomeration makes the size of the characteristic length Tc 
larger than that of non-agglomeration state, it is expected that p^(7agg) > 
^(7agg) fo^ 7agg < 7agg- Here, p^(7agg) represents the statistical average of 
Pc(7agg5'^s)'s as a function of 7agg. In fact. Figure [7] shows that p^(7agg) > 

TchLgg) for 7agg < 7lgg- _ 

When an individual pdlaggyis) becomes smaller than Pc(7agg) for a large 
7agg, the non- vanishing total conductivity crtota.i{p,1a.gg,is) is very small at 
p e (pc(7agg,'^s),Pc(7agg)) and increases weakly with respect to p there. It 
means that the exponent t(7agg,'^s) becomes large in the conductivity curve 
03]). Figure [TTT b) illustrates the correlation between the exponent t(7agg,'^s) 
and the threshold Pc(7agg, is) in ACPMs, which means that the smaller the 
thresholds Pc(7agg, is) are, the larger the exponents t(7agg, is) are. These 
properties are the same as Figure [TIT a). 

However the range of Figure fTTlfb) quite differs from Figure [TlT a). The 
fluctuations in Figure [TO] look symmetric with respect to their averages 
whereas those in Figure [7] are asymmetry. The larger 7agg is, the smaller 
the average of thresholds p^(7agg) is and the larger that of the exponents 
t(7agg) is. Thus the trend of Figure [7| could not be interpreted only by the 
size effect. The difference should be regarded as an shape effect of the ag- 
glomerated clusters. 

In the previous work [1|, we investigated the shape effect on the CPMs 
for spheroid. The thinner the spheroid (of oblate case) is, the smaller we 
have thresholds and the larger we obtain exponents, whereas the thicker the 
spheroid (of prolate case) is, the smaller the threshold and the exponent are. 
As mentioned in Ref. [l|, the shape effect can be also interpreted as the broad 
distribution continuum percolation model (BDCPM). Since the spheroids 
with random orientation can be regarded as a distribution (of probability) 
of the local conductivity, CPMs for a shaped object, e.g., spheroid can be 
interpreted as BDCPM from view point of the probability theory. It means 
that such properties of CPMs of the spheroid can be applied to any shape 
problem in CPM including this case. The thinner the shape of particles (or 
clusters) is, the smaller the threshold is and the larger the exponent is. 

As showed in Subsection I4.1.2[ it could be regarded the larger 7agg is, the 
thinner the shape of the agglomerated cluster becomes. Hence the properties 
of the averaged values in Figure [7] can be interpreted as the shape effects. In 
other words, the tendency of the threshold and the exponent agrees with that 
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of thin shape effects of the spheroids. With the arguments in Subsections 
14.1.11 and I4.1.2[ it means that the larger 7agg is, the smaller we have the 



threshold and the larger we obtain the exponent. 
Figure [7] and Table [TJ 



This properties reproduces 



4-2.3. The universality on conductivity in ACPM 

The fluctuations of the threshold in CPM are essentially the same as those 
of ACPM if the agglomeration parameter 7agg reads the size of p appropri- 
ately. Though the dependence of the average of the exponents on 7agg differs 
from the p dependence, we regard that the difference comes from the fact 
that the agglomerated cluster has a thin shape as another geometrical effect. 

Let us consider the universal properties of ACPMs or the behavior of 
ACPMs in the infinite region. 

Figure [12] exhibits the dependence of the threshold and the exponent on 
the size of the system for the 7agg = 0.5 case. In the computations of Figure 
[T2| we used larger N = N^ = Ny = N^ analyzed regions as A^ = 272 and 
A^ = 344. The relative particle sizes (the radii) are p = 0.79 and p = 0.63. 
Since these computations are harder than A^ = 216, we computed only twenty 
curves for these additional cases respectively as shown in Figure [T^ and Table 
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Figure 12: The dependence of the thresholds and exponents on the size of the system for 
ACPM with 7agg = 0.5. The fiUed circle corresponds to each computation for a seed and 
the dotted line is their average. 



Figure [12] and Table [5] show that the smaller the radius p is, the smaller 
the fluctuations are. This means that the fluctuations comes from the size 
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Tabic 5: The p dependence of the threshold and the exponent for the case 7agg = 0.5. 





Threshold 


Exponent 


p 


Average 


Maximum 


Minimum 


Average 


Maximum 


Minimum 


1.00^ 
0.79* 
0.63* 


0.250 
0.251 
0.253 


0.283 
0.269 
0.266 


0.231 
0.234 
0.243 


1.689 
1.702 

1.714 


1.758 
1.746 
1.744 


1.597 
1.653 
1.676 



*: twenty curves, f thirty curves. 



effect. They also show that the asymptotic values of the thresholds p^ and 
the exponents t differ from those of 7agg = 0.0, and they do not approach to 
those of 7agg = 0.0. We conjecture that the agglomeration, at least in the 
case of our algorithm of agglomeration, has effects on the universality class. 
The difference comes from the shape effect of the agglomerated clusters as 
argued above. 



5. Summary 

By employing the simple algorithm to simulate an agglomerated con- 
figuration of particles, which is controlled by the agglomeration parameter 
7agg ^ [0, 1], we numerically show how the thresholds Pc and the exponent 
t of the conductivity curve ([3]) depend upon the agglomerations as shown 
in Figure [71 The larger the agglomeration parameter 7agg is, the larger the 
fluctuations are. Further the larger 7agg is, the smaller the average of the 
thresholds is and the larger the average of the exponents is. 

From Section 4, we conclude that the origin of these effects could be 
interpreted as the size effect mainly, because the agglomeration makes the 
percolation clusters larger than those of non-agglomerated case. Since for 
the finite region, the enlargement of the clusters makes the fluctuation of the 
conductivity large, the size effect is crucial if the size of system is not large 
enough. 

Subsidiarily the shape of the agglomerated clusters also affect the conduc- 
tivity. It means that the shape of the agglomerated clusters has effects on the 
universality class of the threshold and the exponents as shown in Subsection 
1^:2:31 

Due to the development of the technology, devices with small sizes and 
the small conductive nanoparticles are concerned. Small particles basically 
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agglomerate much due to their surface energy. Then in the conductivity in 
the composite materials of the conductive nanoparticles, the agglomeration 
causes the fluctuations of the conductivity and the difference of individual 
devices. Hence we believe that our results shed a new light on the applications 
of the percolation theory to such a real material system. 
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